Vortex lattice structures of Sr 2 Ru04 
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The vortex lattice structures of Sr2Ru04 for the odd parity representations of the superconducting 
state are examined for the magnetic field along the crystallographic directions. Particular emphasis 
is placed upon the two dimensional representation which is believed to be relevant to this material. 
It is shown that when the zero-field state breaks time reversal symmetry, there must exist two 
superconducting transitions when there is a finite field along a high symmetry direction in the basal 
plane. Also it is shown that a square vortex lattice is expected when the field is along the c-axis. 
The orientation of the square lattice with respect to the underlying ionic lattice yields information 
as to which Ru 4d orbitals are relevant to the superconducting state. 
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^ . The oxide superconductor Sr2Ru04 is structurally similar to the high T c materials but differs markedly from the 
latter in its electronic structure In particular, the normal state near the superconducting transition of Sr2Ru04 is 
O j. well described by a quasi-2D Landau Fermi liquid || . There now exists considerable evidence that the superconducting 
state of Sr2RuC>4 JO is not a conventional s-wave state. NQR measurements show no indication of a Hebel-Slichter 
peak in 1/T\T T c is strongly suppressed by non-magnetic impurities B, and tunneling experiments are inconsistent 
with s-wave pairing ||. While these measurements demonstrate that the superconducting state is non s-wave, they 
do not determine what pairing symmetry actually occurs in this material. The determination of the pairing symmetry 
in unconventional superconductors is a notoriously difficult problem and theoretical insight provides a useful guide. 
The observations that the Fermi liquid corrections due to electron correlations are similar in magnitude to those found 
in superfluid 3 He and that closely related ruthenates are itinerant ferromagnets have led to the proposal that the 
superconducting state in Sr2Ru04 is of odd-parity Even with this insight there still remain five odd-parity states 
that have different symmetry - all of which have a nodeless gap and therefore similar thermodynamic properties ||. 
Recently, /zSR measurements indicate that a spontaneous internal magnetization begins to develop at T c [Q. The 
most natural interpretation of this magnetization is that the superconducting state breaks time reversal symmetry 
i (T). This places a strong constraint on the pairing symmetry in Sr2Ru04 since it implies that the superconducting 
' order parameter must have more than one component ||. Of the possible representations (REPS) of the D±h point 
group only the two dimensional (2D) T$ u and T^g REPS exhibit this property. Of these two the r§ u REP is the most 
likely to occur in Sr 2 Ru04 due to arguments of Ref. [|| and the quasi-2D nature of the electronic properties. The 
order parameter in this case has two components (771, 772) that share the same rotation-inversion symmetry properties 
I as (fcjc, k v ) j|. The broken T state would then correspond to (771, 772) oc (1, i). 

I investigate within Ginzburg Landau (GL) theory the vortex lattice structures expected for the odd-parity REPS 
of the superconducting state; focusing mainly on the F5„ REP. It is initialy shown that a general consequence of 
the broken T state described above is that in a finite magnetic field oriented along a high symmetry direction in 
the basal plane there will exist a second superconducting transition in the mixed state as temperature is reduced. 
The observation of such a transition would provide very strong evidence that the order parameter belongs to the T§ u 
REP. The form of the vortex lattice for a field along the c-axis is then investigated for both the one dimensional 
(ID) and the 2D T^u REPS of the superconducting state. It is shown that a square vortex lattice is expected to 
appear for all the REPS, however observable differences exist between the ID and the 2D REPS. Finally, within the 
recently proposed model of orbital dependent superconductivity of Sr2Ru04 it is also shown that the orientation 
of square vortex lattice with respect to the underlying crystal lattice dictates which of the Ru 4d orbitals give rise to 
the superconducting state. 

To demonstrate the presence of the two superconducting transitions described above consider the magnetic field 
along the x direction (x is chosen to be along the basal plane main crystal axis) and a homogeneous zero-field state 
(77!, 772) cx (1, i). In general the presence of a magnetic field along the x direction breaks the degeneracy of the 
(^1,^2) components, so that only one component will order at the upper critical field [e.g. (771,7/2) oc (0,1)]. As 
has been shown for type II superconductors with a single component the vortex lattice is hexagonal at T c and the 
order parameter solution is independent of x so that o~ x (a reflection about the x direction) is a symmetry operation 
of the (771,772) cx (0,1) vortex phase flipf . Now consider the zero-field phase (771,772) cx cr x transforms (l,i) to 

(—l,i) 7^ e l< ^(l,i) where 4> is phase factor. This implies that o~ x is not a symmetry operator of the zero-field phase. It 
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follows that there must exist a transition in the finite field phase at which 771 becomes non-zero. Similar arguments 
hold for the field along any of the other three crystallographic directions in the basal plane. Evidence for this transition 
may already exist in the ac magnetic susceptibility measurements of Yoshida et. al. They observed a second 

peak in the imaginary part of the magnetic susceptibility only when the flux lines were parallel to the basal plane. 
For a more detailed analysis consider the following dimensionless GL free energy density for the T^ u REP 

/ = - M 2 + hT/2 + Mmv* 2 - mvl) 2 /2 + /%M 2 M 2 + \d xVi \ 2 + \D yV2 \ 2 (1) 

+ K 2 (\D yVl \ 2 + \D xV2 \ 2 ) + k 5 (\D zVi \ 2 + \D z m?) 
+ K 3 [{D xm )(D yV2 )* + h.c] + Ki[(D ym ){D xV2 )* + h.c] + h 2 . 

where h = V x A, D v = V„/k — iA u , f is in units B 2 /(Att), lengths are in units A = [h c 2 0\j (16e 2 Kian)] 1 ^ 2 , h 
is in units y/2B c = $o/(4ttAC) 1/2 , a = a (T - T c ), £ = (ni/a) 1/2 , and k = A/£. Note that A, £, B c and k are 
simply convenient choices and do not correspond to measured values of these parameters. A thorough analysis of 
Eq. |l| is difficult due to the unknown phenomenological parameters (3 2l (3^, k 2 , K3, and K4. To simplify the analysis 
these parameters are determined in the weak-coupling, clean-limit for an arbitrary Fermi surface. Taking for the 
REP the gap function described by the pseudo-spin-pairing gap matrix: A = i\j]\v x j y/ (v 2 ) + rj 2 v y / ^/ (v x )]<J z (Ty 
where the brackets () denote an average over the Fermi surface and oi are the Pauli matrices, it is found that 
P2 = K 2 = = K4 = 7 and 0a = 37 — 1 where 7 = (v x v 2 )/(v x ). Note that < 7 < 1 and that 7=1/3 for a 
cylindrical or spherical Fermi surface. These parameters agree with the cylindrical Fermi surface results of Ref. |fl2| . 
It is easy to verify that in zero-field (771, 772) <x (1, i) is the stable ground state for all 7. 

It is informative to determine the values of 7 that are relevant to Sr 2 Ru04. LDA band structure calculations 
p3| , p^| reveal that the density of states near the Fermi surface are due mainly to the four Ru 4cZ electrons in the 
t 2g orbitals. There is a strong hybridization of these orbitals with the O 2p orbitals giving rise to antibonding ir* 
bands. The resulting bands have three quasi-2D Fermi surface sheets labeled a, (3, and 7 (see Ref. pi). The a and (5 
sheets consist of {xz,yz} Wannier functions and the 7 sheet of xy Wannier functions. In general 7 is not given by a 
simple average over all the sheets of the Fermi surface. A knowledge of the pair scattering amplitude on each sheet 
and between the sheets is required to determine 7 |)|,[l5| . Recently, to account for the large residual density of states 
observed in the superconducting state, it has been proposed that either the xy or the {xz,yz} Wannier functions 
exhibit superconducting order This model implies that that there are two possible values of 7; one for the 7 
sheet ("fxy) and one for an average over the {a, j3} sheets (jxz.yz)- A tight binding model indicates j xy = 0.67 and 
lxz,yz = 0.11 |l6| |. These values are sensitive to changes in the parameters of the tight binding model, however the 
qualitative result that "f xy > 1/3 and j xz yz < 1/3 is robust. Physically "f xy > 1/3 because of the proximity of the 7 
Fermi surface sheet to a Van Hove singularity and "/ X z, y z < 1/3 due to quasi ID nature of the {a,/3} surfaces p3| , pT[ . 

Following Burlachov jlTj for the solution of upper critical field for the field in the basal plane, the vector 
potential is taken to be A = Hz(s'm 8, — cos 9, 0) {9 is the angle the applied magnetic field makes with the x direction). 
After setting the component of D along the field to be zero it is found that H%{0) = k/(k 5 A(6»)/2) 1 / 2 where A(0) = 
1 + 7 — [(1 — 7) 2 — (1 + 7)(1 — 37) sin 2 29] 1 / 2 . A measurement of the temperature independent four-fold anisotropy 
in thus determines 7. To determine the field at which the second transition discussed above occurs consider the 
magnetic field along the x direction. The free energy of Eq. [j] is then similar to that studied in UPt3 p^Jl9| , p^| and 
since Si'2Ru04 is a strong type II superconductor with a GL parameter of 31 for the field in the basal plane T\Jm the 
procedure of Garg and Chen |l8) to study the second transition can be applied here. At r\\ orders and and the 
vortex lattice solution is given by [|l0|Jl^Jl9|l 

T) 1 =y2 Cr[ e^ e -(' t E/7) 1/2 «Hfe-«n/(«H)] 2 /2_ ( 2 ) 

where c n = e m27r / 2 and q has the two possible values qf = \/3H ktt^ / K5) 1 / 2 or q 2 — H ktt(j / K5) 1 / 2 / V3 (these two 
solutions are degenerate in energy). At the second transition the rj 2 component becomes non-zero. As is discussed 
in Refs. Jl^Jl^] the solution for r\ 2 corresponds to a lattice that is displaced relative to that of r\\ by d = (y,z). 
Accordingly, the field at which the second transition occurs is found by substituting 

T\ 2 = ir ^ Cne i(ng+KHy)(z-z) e -^ r J^KH[y-y-qn/(KH) 2 ]/2 ^ 

and Eq. [2] into the free energy, minimizing with respect to the displacement vector d, and determining when the 
coefficient of r 2 becomes zero. This yields for the ratio of the second transition (H 2 ) to the upper critical field 
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H 2 _ 1/2 (3 A - 7(25! - \S 2 \) m 
Hg 7 /?a-7 3/2 (2Si-|S2|) 



where Pa = 1-1596, Si = |?/i| 2 |?72| 2 /(|?7i| 2 |?72| 2 ), S 2 = ( ? ?i J /2) 2 /(l 7 ?i| 2 l^l 2 ), the over-bar denotes a spatial average, 
and the subscript min means take the minimum value with respect to d and with respect to q — qi or q — q 2 . The 
numerical solution of Eq. |4| is shown in Fig. [j]. Three vortex lattice configurations are found to be stable as a function 
of 7 (depicted in Fig. |l]). For < 7 < 0.187 q = q 2 and d = (T y ,T z )/A (T y and T z are the translation vectors of the 
centered rectangular cell for the r\\ lattice), for 0.187 < 7 < 0.433 q = qi and d = (T y ,T z )/4, and for 0.433 < 7 < 1 
q = q 2 and d = 0. For the field along x ± y the ratio H 2 /H C2 is given by Eq. || with 7 replaced by (1 — 7)/(l + 37). 
As a consequence the form of the vortex lattice will depend on the field direction. As has been discussed in detail in 
Ref. [pTJ the shape of the vortex lattice unit cell for H < H 2 will be strongly field dependent. 
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FIG. 1. The ratio of the two transition lines for the field along x as a function of 7. The open (closed) circles correspond 
to the zeroes of the 772 (j?i) lattice. The vertical lines separate regions where the depicted vortex lattice structures are favored. 
For all three lattice structures the y and z axes have the same orientation and the dimensions of the rectangular cell are the 
same. 

Consider the magnetic field oriented along the c-axis. Setting D z = writing II + = k(iD x + D y )/(2H), II_ = 
k(iD x — D y )/(2H), rj + = (r) x + ii] y )/\/2, and r?_ = {rj x — ir] y )/\/2, minimizing the quadratic portion of Eq. [I] with 
respect to 77+ and r\- yields 

oJ^A-rrf (l+7)(l + 2iV) (i + 7)n 2 + + (i-3 7 )n 2 \ f v+ \ 

UJ V(! + 7)n 2 _ + (1-37)14 (1+ 7 )(1 + 2A) )\V-J- [) 

where N = n + n_. The maximum value of H that allows a non-zero solution for (rj +1 r)_) yields the upper critical 
field H^ 2 ■ For 7 ^ 1/3 H£ must be found numerically (note that for 7 = 1 /3 the solution can be found analytically 
pl|,p|). Expanding (r/ + ,r;_) in terms of the eigenstates of N (Landau levels) up to N = 32 and diagonalizing the 
resulting matrix yields the result for (7) shown in Fig. |2[ The solution for the form of the vortex lattice represents 
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a complex problem due to presence of many Landau levels in the solution of (rj + , ) and the weak type II nature 
of Sr2Ru04 for the field along the c axis (Ref. indicates k ~ 1.2). Here I present results that are strictly valid 
in the large k limit and leave the treatment for general k to a later publication (a perturbative calculation indicates 
that the qualitative results are unchanged for k — 1.2) |l6). The form of the eigenstate of the H£ solution is 
found to be rj + (r) = J]„>o a 4n+2?W+2(r) and r)-(r) = Y, n >o «4n^4n(r) where 0„(r) = Y, m c m e lqmv 2- n / 2 H n {x - 
qm/(KH))e~ KH ( x ~ qm /( KH ^ / 2 /{n\) 1 / 2 where the coefficients a n are real, (x, y) is the vector (x,y) rotated by an angle 
9 about the z axis and H n (x) represent Hermite polynomials. In the large k limit the form of the vortex lattice is 
found by minimizing 

/3 = h/(M 2 ) 2 (fi is the quartic part of Eq. |) with respect to the coefficients c„, q, and 9. It 
is assumed that c„ = c„ +2 - This restricts the vortex lattice structures to be centered rectangular with a short axis 
L y = 2ir/q and a long axis L x = 2q/(nH). The ratio t = L x /L y is \/3 for a hexagonal vortex lattice and is 1 for 
a square vortex lattice. I further restrict the analysis to the two orientations 9 = {0, 7r/4} since these correspond to 
aligning one of the vortex lattice axes with one of the high symmetry directions in the basal plane. Remarkably, the 
treatment of the many Landau levels in the solution of rj+ and rj- becomes numerically straightforward when (3 is 
expressed as a sum over the reciprocal lattice given by 1 = xl\2ir/L x + y^n/Ly (see also Ref. Q). It is found 
that (3 is minimized for c„ = e m ~*l 2 and that the values of t and 9 depend upon 7. For 7 < 1/3 (7 > 1/3) 9 = 
(9 = 7r/4) and t varies continuously from \/3 to 1 as 7 decreases (increases) from 1/3 to 1/3 — 0.0050 (1/3 + 0.0050). 
For 7 < 1/3 — 0.0050 and 7 > 1/3 + 0.0050 the minimum (3 corresponds to t = 1. This implies that for "f X z, V z a square 
vortex lattice rotated 7r/4 about the c axis from the crystal lattice is expected and for 7^ a square vortex lattice that 
is aligned with the underlying crystal lattice is expected near H£ 2 ■ Note the appearance of the square vortex lattice 
correlates with an anisotropy in of |1 - H^(9 = 0)/H^(9 = tt/4)| > 0.01. 
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FIG. 2. The inverse of the upper critical field as a function of 7 for the field applied along the c-axis. 

The recent observation of a square vortex lattice in Sr2Ru04 |p3[ makes it of interest to compare the above behavior 
to that expected for the ID REPS of D^. It is well known that for single component order parameters non-local 
corrections to the standard GL theory stabilize a square vortex lattice |24|-|26|] . In particular the the following non-local 
term will stabilize the square lattice 

e[\(D 2 x ~ D 2 y )^\ 2 - \(D x D y + D y D x )^\ 2 }. (6) 

Treating this term as a perturbation to the GL free energy leads to ijj — 4>o — £04 where e = \/oeH/K (k is the GL 
parameter) near H£ . As e increases (note e = at T c ) the vortex lattice continuously distorts from hexagonal to 
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square |24|-|26| until the square vortex lattice is stable for |e| > 0.024. The sign of e determines the orientation of the 
vortex lattice; for e > the vortex lattice is aligned with the underlying lattice while for e < the lattice is rotated 7r/4 
with respect to the underlying crystal lattice 24 2q|. The sign of e has been determined within a weak coupling clean 
limit approximation for the ID odd-parity REPS. For the Ai u REP using A = ip(xv x / y/ (v%) + yv y / y/ (i>|)) ' Sia v the 
sign of e is determined by the sign of 3((w^ + Vy)v x Vy) j '((v x + Vy)v x ) — 1. Using a form for A that is analogous to that 
used for the A\ u REP, the same result is found for all the ID odd-parity REPS. This implies that the final orientation 
of the square vortex lattice for the ID REPS is the same as that found for the 2D REP for superconducting order in 
the xy or the {xz,yz} orbitals. The behavior of the vortex lattice for the ID REPS as a function of e is very similar 
to that for the 2D REP as a function of 7. An observable difference between the 2D and the ID REPS is that for 
the 2D REP the vortex lattice remains square up to T c while for the ID REPS the vortex lattice is hexagonal at T c . 
Also, the GL theories for the ID and the 2D REPS predict a four-fold anisotropy in but this anisotropy vanishes 
at T c for the ID REPS and does not vanish at T c for the 2D REP. 

In conclusion I have examined GL models for the odd-parity REPS of the superconducting state for Sr2Ru04. It 
was found that if the zero-field ground state breaks T symmetry (the 2D REP) then there should exist a second 
transition in the mixed state when the magnetic field is applied along a high symmetry direction in the basal plane. 
It was also shown that when the field is along the c-axis there will be a square vortex lattice for all the possible odd 
parity superconducting states. 
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